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Abstract — During  the  2002,  2003,  and  2004  Hurricane  seasons,  and 
the  2005  Ocean  Winds  and  Rain  Winter  Experiment,  the  University  of 
Massachusetts  (UMass)  installed  two  instruments  on  the  NOAA  N42RF 
WP-3D  research  aircraft:  the  Imaging  Wind  and  Rain  Airborne  Pro¬ 
filer  (IWRAP)  and  the  Simultaneous  Frequency  Microwave  Radiometer 
(SFMR).  IWRAP  is  a  dual-band  (C-  and  Ku),  dual-polarized  pencil- 
beam  airborne  radar  that  profiles  the  volume  backscatter  and  Doppler 
velocity  from  rain  and  that  also  measures  the  ocean  backscatter  response. 
It  simultaneously  profiles  along  four  separate  incidence  angles  while 
conically  scanning  at  60  RPM.  SFMR  is  a  C-band  nadir  viewing 
radiometer  that  measures  the  emission  from  the  ocean  surface  and 
intervening  atmosphere  simultaneously  at  six  frequencies.  It  is  designed 
to  obtain  the  surface  wind  speed  and  the  column  average  rain  rate.  Both 
instruments  have  previously  been  flown  during  hurricane  seasons  2002, 
2003  and  2004. 

Spectral  techniques  can  been  used  to  differentiate  the  contributions 
of  the  volume  backscatter  from  rain  from  the  ocean  surface  backscatter. 
Initial  results  from  the  Winter  2005  datasets  are  presented.  Attenuation 
from  rain  is  also  particularly  important  at  high  microwave  frequencies 
such  as  Ku-band.  Thanks  to  IWRAP’s  profiling  capabilities,  Ku-band 
attenuation  models  and  dropsize  distribution  (DSD)  parameters  of  the 
observed  precipitation  from  dual-wavelength  techniques  within  tropical 
cyclones  are  presented. 

I.  Introduction 

IWRAP,  the  Imaging  Wind  and  Rain  Airborne  Profiler,  is  a  high- 
resolution,  dual-band  (C-  and  Ku-band),  dual-polarized,  pencil-beam 
radar  that  profiles  the  Doppler  velocity  and  volume  backscatter  from 
rain  and  measures  the  ocean  backscatter  response  at  15  to  120  m 
range  resolution.  [1].  Developed  by  the  University  of  Massachusetts 
with  the  support  of  Remote  Sensing  Solutions,  and  funded  by 
NOAA/NESDIS,  NASA/JPL  and  ONR,  the  system  operates  simul¬ 
taneously  at  four  separate  incidence  angles  (covering  an  effective 
continuous  swath  of  25  to  50  degree  incidence)  while  conically 
scanning  at  60  RPM.  SFMR,  the  Simultaneous  Frequency  Microwave 
Radiometer,  is  a  C-band  nadir  vieweing  radiometer  that  measures  the 
emission  from  the  ocean  surface  and  atmosphere  simultaneously  at 
six  separate  frequencies  (approximately  ranging  from  4  to  7  GHz).  It 
is  used  operationally  to  provide  continuous  estimates  of  the  surface 
wind  speed  and  the  column  average  rain  rate  [2]. 

During  the  2001,  2002  and  2003,  the  UMass  IWRAP  instrument 
was  respectively  built,  tested,  and  deployed  aboard  the  NOAA 
WP-3D  hurricane  research  and  reconnaissance  aircraft.  Since  then, 
IWRAP  has  participated  every  year  in  the  NOAA/NESDIS  Hurricane 
Ocean  Winds  and  Rain  Experiments,  the  NOAA/NESDIS  Winter 
Experiments,  and  the  ONR  CBLAST  Hurricane  Program.  During 
these  fields  experiments,  dozens  of  missions  have  been  flown  through 
winter  extratropical  storms  at  high  latitudes,  as  well  as  tropical  storms 


and  hurricanes  (such  as  Gustav,  Isidore  and  Lili  in  2002,  Fabian  and 
Isabel  in  2003,  and  Frances,  Ivan  and  Jeanne  in  2004),  during  various 
phases  of  these  systems.  All  the  IWRAP  data  was  acquired  with 
coincident  SFMR  measurements  and  hundreds  of  dropsonde  and  buoy 
in- situ  measurements.  In  2005,  the  IWRAP  instrument  was  upgraded 
to  acquire  raw  data. 

In  this  paper,  we  present  initial  results  from  these  datasets,  which 
establish  the  grounds  to  ultimately  analyze  and  quantify  the  impact 
of  rain  on  the  ocean  surface  backscatter  response.  The  structure 
is  as  follows:  we  first  present  the  use  of  spectral  techniques  with 
the  2005  Winter  data,  which  allows  to  individuate  and  separate  the 
contributions  of  the  volume  backscatter  from  rain  from  the  ocean 
surface  backscatter.  The  second  section  describes  the  use  of  the 
coincident  dual-wavelength  measurements  to  derive  an  attenuation 
model  and  rainfall  rate  estimates.  Finally,  the  last  section  uses  these 
results  to  derive  a  two-parameter  drop- size  distribution  within  tropical 
cyclones. 

II.  Ocean  and  rain  spectra 

The  IWRAP  system  collected  raw  data  for  the  first  time  during  the 
2005  winter  experiment.  The  purpose  of  the  upgrade  was  to  address 
one  of  the  fundamental  limitations  of  pulse-pair  processing,  which 
was  implemented  on  the  FPGAs  of  the  instrument’s  data  acquisition 
system.  The  limitation  comes  from  the  incorrect  estimation  of  the 
Doppler  moments  (mean  and  spectral  width)  when  more  than  one 
contribution  exists  in  the  Doppler  spectrum.  For  the  case  of  any 
radar  profiling  the  atmosphere  and  the  ocean  surface,  the  responses 
from  volume  backscatter  from  rain,  surface  backscatter  from  the 
ocean  surface,  and  nadir  returns  (in  non  nadir-looking  systems)  can 
simultaneously  contribute  to  the  measurements  [3]  [4].  Given  the 
fact  that  such  contributions  will  typically  be  centered  in  different 
positions  within  the  span  of  the  full  Doppler  spectrum,  the  pulse-pair 
estimates  will  estimate  a  Doppler  velocity  mean  which  is  a  power- 
weighted  combination  of  the  individual  Doppler  velocity  means  of 
each  one  of  the  contributions.  The  spectral  measurements  provide 
instead  the  means  to  individuate  the  different  contributions.  Moreover, 
the  ocean  spectral  response  can  be  modelled  under  all  wind  conditions 
and  in  both  rain  and  rain-free  conditions.  Rain  processes  can  be 
characterized  from  the  dual- wavelength  measurements,  as  it  will 
be  discussed  in  the  next  section.  The  exact  knowledge  of  both 
contributions  thus  allows  to  precisely  quantify  the  impact  of  rain 
on  the  ocean  backscatter  on  the  spectral  response  as  well  as  in  the 
integrated  power. 
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Fig.  1.  Doppler  spectra  at  altitudes  ranging  from  500  to  30  m  at  a  fixed 
azimuth  angle  are  shown.  The  different  contributions  from  rain,  nadir  return 
and  ocean  surface  backscatter  are  clearly  indicated. 


Fig.  2.  Doppler  spectra  for  half  an  azimuthal  scan  at  two  altitudes  (500  and 
30  m)  are  shown.  The  different  contributions  from  rain  and  ocean  surface 
backscatter  are  indicated. 


Figure  1  shows  the  derived  Doppler  spectra  at  altitudes  ranging 
from  30  to  500  m  at  a  fixed  azimuth  angle  (corresponding  to 
the  downwind  look  direction)  and  one  incidence  angle  (40  degrees 
incidence)  in  the  presence  of  light  rain  (less  than  10  mm  hr-1) 
and  wind  speeds  below  20  m-s_1.  The  figure  clearly  shows  that 
the  rain  contribution,  shown  alone  at  an  altitude  of  510  m,  is  first 
contaminated  by  the  nadir  return  at  210  m,  and  then  by  the  ocean 
surface  backscatter  below  120  m.  In  this  situation  it  is  clear  that 
the  nadir  contributions  can  be  easily  removed  by  filtering  out  the 
undesired  part  of  the  spectrum  for  the  lowest  boundary  layer.  The 
relatively  low  winds  place  both  ocean  and  rain  spectral  components 
close  enough  as  to  make  the  filtering  a  non  viable  solution.  However, 
the  ocean  response  at  0  m  above  the  ocean’s  surface  can  be  used  to 
model  the  responses  at  other  altitudes  to  substract  them  from  the 
obtained  spectrum.  This  approach  is  currently  being  evaluated.  In 
the  presence  of  stronger  winds,  like  the  ones  found  in  the  context  of 
tropical  cyclones,  both  contributions  are  expected  to  be  placed  further 
apart  from  each  other,  and  thus  filtering  could  be  used,  although  only 
for  some  azimuthal  positions  along  the  scan;  as  figure  2  shows,  the 
ocean  backscatter  response  (shown  at  30  m,  where  it  predominates 
over  the  contribution  from  the  rain)  overlays  the  rain  spectra  for  some 
of  the  azimuthal  positions.  Depending  on  the  technique  employed, 
however,  the  rain  spectra  does  not  need  to  be  individuated  at  all 
azimuthal  positions  in  order  to  retrieve  the  boundary  layer  winds. 

III.  DUAL- WAVELENGTH  ATTENUATION  FROM  PRECIPITATION 

Many  remote  sensing  techniques  have  been  developed  for  precipi¬ 
tation  rate  ( Rr )  estimation.  These  methods  can  be  broken  down  into 
single  and  multi-parameter  approaches.  One  of  the  most  commonly 
used  single  parameter  techniques  relates  the  radar  reflectivity  factor  Z 
(i.e.  volume  backscatter  measurements)  to  Rr,  usually  by  means  of  a 
power  law.  Unfortunately,  the  Z-Rr  relationship  can  vary  significantly 


for  different  types  of  precipitation  and  for  different  climates.  Battan 
lists  more  than  60  different  Z-Rr  relationships  based  on  drop-size 
distribution  measurements  made  in  different  climatic  regions  [9].  To 
assume  that  is  it  possible  to  find  such  a  Z-R  relationship  is  to  imply 
that  some  of  the  drop- size  distribution  (DSD)  parameters  are  not 
independent  variables. 

Multi-parameter  approaches  are  better  suited  to  retrieve  precip¬ 
itation  rate.  One  approach  is  to  acquire  volume  backscatter  mea¬ 
surements  at  two  different  electromagnetic  wavelengths,  a  long 
wavelength  (in  our  case,  C-Band)  that  is  virtually  not  attenuated, 
and  a  shorter  wavelength  (Ku-band)  that  is  attenuated.  From  both 
these  measurements,  the  path  averaged  specific  attenuation,  Ks ,  in 
dB  km-1,  can  be  determined  from  the  reflectivity  factors  at  the  two 
frequencies  and  over  two  different  ranges,  as  given  by: 

Ks  ~  AT - v  [Zc(t2)  +  Zku(ri)  —  (Z(Ti)  —  Zku{r 2))]  .  (1) 

2  [r2  —r  1)  L 

where  n  and  r 2  are  the  ranges  and  Zc(ri)  and  Zku{ri)  denote 
the  reflectivity  factors  in  dBZ  at  range  ri  at  C-  and  Ku-band, 
respectively.  From  measurements  of  Ks  and  Z,  a  two  parameter 
drop-size  distribution,  N(D ),  can  be  calculated,  and  from  this 
information,  the  precipitation  rate  can  be  determined.  Note  that  such 
a  method  still  relies  on  measuring  Z,  but  provides  another  degree 
of  freedom  to  better  describe  the  drop- size  distribution  that  is  being 
observed.  Moreover,  the  attenuation  estimate  is  independent  of  the 
radar  calibration.  However,  it  requires  that  no  hail  is  present  between 
r  1  and  r2\  otherwise,  Ks  will  remain  the  average  attenuation  along 
the  path. 

Figure  3  presents  a  scatter  plot  of  the  specific  attenuation  versus 
the  SFMR  rain  rate  estimates.  The  differential  reflectivity  is  obtained 
from  reflectivity  measurements  at  flight  level  and  1  km  range  distance. 
Two  models  derived  by  Atlas  [5]  at  10°C  and  40°C  are  overlaid.  From 
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Fig.  3.  Scatter  plot  of  specific  attenuation  measurements  from  dual¬ 
wavelength  techniques  obtained  during  three  flights  through  hurricane  Isabel 
(2003)  versus  the  SFMR  rainfall  rate  estimates.  Atlas’  models  of  the  specific 
attenuation  at  10°C  and  40° C  are  overlaid.  The  derivation  of  the  corrected 
SFMR  rain  rate  is  described  in  [8]. 


these  measurements,  a  power  law  relating  the  specific  attenuation  and 
the  SFMR  rate  can  be  derived.  The  obtained  relationship  is: 

Ks  =  0.02i^'276.  (2) 

The  1-cr  errors  for  the  fitted  parameters  are  lxlO-3  for  the  coefficient 
and  1.4xl0-2  for  the  exponent.  The  derived  K-R  relationship  can 
now  be  used  to  obtain  rainfall  rates  from  the  dual- wavelength  IWRAP 
reflectivity  measurements.  Figure  4  presents  an  instance  of  such 
measurements  corresponding  to  a  flight  through  Hurricane  Isabel 
(2003).  As  the  figure  shows,  the  agreement  is  excellent  except  at 
the  highest  rain  rates  (>  70  mm  hr-1),  where  the  IWRAP  estimates 
predict  a  greater  rainfall  rate  than  SFMR.  Possible  reasons  for  the 
discrepancy  are  associated  to:  1)  the  SFMR  rainfall  rate  estimates  are 
column  averaged,  while  the  IWRAP  uses  only  the  attenuation  within 
the  first  km  from  the  aircraft  in  deriving  the  estimates.  These  volumes 
are  different  and  for  the  most  intense  rain  rates  this  difference  can 
be  significant.  Its  cone  geometry  makes  SFMR  more  sensitive  to  rain 
rates  closer  to  the  surface;  and  2)  IWRAP  is  a  much  higher  resolution 
instrument,  and  it  is  likely  to  see  higher  rain  rates  since  SFMR  will 
smooth  the  precipitation  events  and  miss  high  resolution  events. 

IV.  Two-parameter  Drop-size  distributions 

By  using  both  the  attenuation  and  reflectivity  measurements,  a 
two-parameter  DSD  can  be  derived.  We  will  hereafter  consider  an 
exponential  distribution  of  the  form 

N(D)  =  N0e~AD,  (3) 

where  N0  and  A  are  functions  of  the  diameter  of  the  hydrometers, 
D.  With  this  assumption,  the  specific  attenuation  can  be  expressed 
by 

Ka  =  4.34*103A^o  [  e~ADae(D)dD,  (4) 

J  D 


where  Ks  is  in  dB-km-1  and  ae  is  the  total  attenuation  or  extinction 
cross  section.  The  expression  can  be  approximated  by  a  power 
law  [6].  On  the  other  hand,  when  D  is  small  compared  to  the 
electromagnetic  wavelength  A  (i.e.  Rayleigh  scattering,  D  <  A/8), 
the  radar  reflectivity  factor  can  be  expressed  as 

roo 

Z=  N(D)D6dD.  (5) 

Jo 

For  both  expressions,  the  integration  can  be  performed  using  the 
exponential  distribution  given  by  eq.  3.  Figures  5  and  6  show 
the  resulting  A  and  No  parameters  from  the  observations  acquired 
during  three  flights  through  Hurricane  Isabel,  where  the  reflectivity 
measurements  at  C-band  and  flight  level  have  been  used  for  eq.  5. 

iVo  values  have  been  suggested  as  2.2xl04  m_3mm_1  for  rain, 
105  m_3-mm_1  for  snow  and  4xl03-m-3  mm-1  for  graupel  [7]. 
Figure  6  indicates  that  all  these  three  forms  are  present  for  rain  rates 
mostly  below  40  mm-hr_1.  Above  that,  the  median  N0  is  of  the  order 
of  2.8xl04  m_3mm_1,  indicating  the  presence  of  mostly  rain  only. 
The  median  mass  diameter,  Do,  can  then  be  obtained  by  applying 
the  exponential  DSD  under  consideration  into  its  defining  equation: 

[  °  D3  e~AD  dD  =  |  f  D3e~ADdD.  (6) 

Jo  2  J o 

Figure  7  shows  the  derived  median  mass  diameter  for  the  same 
observations. 

These  results  suggest  that  correcting  the  ocean  backscatter  response 
from  attenuation  becomes  an  extremely  difficult  task  even  at  the  lower 
rainfall  rates  (above  10  mm  hr-1),  given  the  large  range  of  raindrop 
diameters.  Even  in  the  situations  where  the  rainfall  rate  is  low  enough 
so  that  the  impact  of  the  attenuation  is  not  large  and  the  impact  of  the 
rain  drops  on  the  ocean  surface  is  believed  not  to  have  too  much  of 
an  effect  for  the  retrieval  of  the  ocean  wind  vector,  the  broad  spread 
in  the  median  diameters  would  make  any  correction  method  prone  to 
large  errors. 

V.  Summary 

For  the  first  time,  IWRAP  raw  data  measurements  of  the  atmo¬ 
spheric  boundary  layer  and  the  ocean  surface  have  been  acquired. 
The  use  of  spectral  techniques  to  identify  and  separate  ocean  and 
rain  backscatters  allows  us  to:  1)  derive  the  ocean  spectral  response 
in  saturation;  2)  analyze  the  impact  of  rain  on  the  ocean  backscatter 
with  an  exact  knowledge  of  both  contributions;  and  3)  remove  the 
ocean  backscatter  response  to  derive  the  ABL  winds  down  to  the 
surface.  Also,  from  the  dual-wavelength  measurements  of  rain,  an 
attenuation  versus  rain  rate  model  has  been  derived,  and  rainfall  rate 


(a) 

Fig.  4.  Derived  rainfall  rate  estimates  from  IWRAP  measurements.  The 
SFMR  estimates  are  also  shown.  The  measurements  correspond  to  an  eye- 
wall  penetration  through  Hurricane  Isabel  on  September  12th,  2003. 
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Fig.  5.  Scatter  plot  of  the  slope  of  the  DSD,  A,  versus  the  corrected  SFMR 
rainfall  rate  estimates  obtained  from  the  C-band  30  degrees  incidence  flight 
level  reflectivity  measurements  acquired  during  the  Isabel  flights  on  the  12th, 
13th  and  14th  of  September. 


estimates  through  hurricane  passes  have  been  presented.  Drop  size 
distribution  parameters  have  also  been  extracted,  indicating  the  diffi¬ 
culties  associated  with  the  correction  of  the  ocean  surface  backscatter 
response  in  the  presence  of  rainfall  rates  above  10  mm  hr-1. 
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Fig.  6.  Scatter  plot  of  the  TVo  coefficient,  No,  versus  the  corrected  SFMR 
rainfall  rate  estimates  obtained  from  the  C-band  30  degrees  incidence  flight 
level  reflectivity  measurements  acquired  during  the  Isabel  flights  on  the  12th, 
13th  and  14th  of  September. 
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Fig.  7.  Scatter  plot  of  the  DSD  slope,  Do,  versus  the  corrected  SFMR 
rainfall  rate  estimates  obtained  from  the  C-band  30  degrees  incidence  flight 
level  reflectivity  measurements  acquired  during  the  Isabel  flights  on  the  12th, 
13th  and  14th  of  September. 
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